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1. Introduction 

Twisted mass fermions [|I|] provide an attractive formulation of lattice QCD that allows for 
automatic &{a) improvement, infrared regularization of small eigenvalues and fast dynamical sim- 
ulations. A particularly attractive feature for the calculation of the nucleon form factors discussed 
in this work is the automatic ^ (a) improvement obtained by tuning only one parameter, requiring 
no further improvements on the operator level. Important physical results are emerging using gauge 
configurations generated with two degenerate flavors of twisted quarks (A^^ = 2) in both the me- 
son and baryon [^] sectors. An example is the accurate determination, using precise results in 
the meson sector, of low energy constants of great relevance to phenomenology. Currently, Nf = 2 
simulations are available for pion mass in the range of about 260-470 MeV for three lattice spacings 
a < 0.1 fm allowing for continuum and chiral extrapolations. In this work we discuss high-statistics 
results on the nucleon form factors obtained at one value of the lattice spacing. Electromagnetic 
and axial form factors (FFs) of the proton and the neutron are fundamental quantities that yield 
information on their internal structure such as their size, magnetization and axial charge. They 
have been studied experimentally for a long time improving their measuements both in terms of 
precision as well as in terms of probing larger momentum transfers. Several lattice collaborations 
are currently using dynamical fermions to calculate these fundamental quantities |5|]. 

The action for two degenerate flavors of quarks in twisted mass QCD is given by 

X 

where we use the tiee-level Symanzik improved gauge action Sg. The quark fields % are in the 
so-called "twisted basis" obtained from the "physical basis" at maximal twist by the transforma- 
tion i/A = -^[1 + ix^'Ys\x and i// = ^-^[1 + /t^/j]. We note that, in the continuum, this action is 
equivalent to QCD. A crucial advantage is the fact that by tuning a single parameter, namely the 
bare untwisted quark mass to its critical value mcr, physical observables are automatically 0'{a) 
improved. A disadvantage is the explicit flavor symmetry breaking. In a recent paper we have 
checked that this breaking is small for baryon observables for the lattice spacing discussed here [^]. 

To extract the nucleon FFs we need to evaluate the nucleon matrix 
elements {N{pf,Sf)\j^\N{pi,Si)), where \N{pf,Sf)), \N{pi,Si)) are 
nucleon states with final (initial) momentum Pf{pi) and spin Sf{si) 
and 7^ is either the electromagnetic current V(f^(x) = ^u{x)y^u{x) — 
\d{x)'Yp,d{x) or the axial current (jc) = !//(;<;) 7^ 75 Whereas 
the matrix element of the axial current receives contributions only 
from the connected diagram shown in Fig. 1 the electromagnetic 
one has, in addition, disconnected contributions. In the isospin limit 
the matrix element of the isovector electromagnetic current V^' [x) = Fig. 1 : Connected nucleon 
W{^)y^'^W{^) has no disconnected contributions [0]. Therefore in three-point function, 
this work we only evaluate the isovector nucleon FFs obtained from 
the connected diagram. 

The electromagnetic matrix element of the nucleon can be expressed in terms of the Dirac and 
Pauli form factors, and F2 defined, in Euclidean time, as 



(1.1) 
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{N{pf,Sf)\V^mN{pt,Si))-. 
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u{pf,Sf)ff^u{pi,Si) 
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N 



EN{Pf)EN{Pi) 

with q = Pf — Pi the momentum transfer and = —q^. These are related to the Sachs electric Ge 
and magnetic Gm FFs via: GsiQ^) =F,{Q^) - ^^^2(2') and Gm{Q^) = F,{Q^) +^2(2')- 

Similarly, the axial current matrix element of the nucleon {N{pf,Sf)\A''^{0)\N{pi,Si)) can be 
expressed in terms of the form factors Ga and Gp with given by 



.2x , ..^''75^ ,^2x 



2" 



2. Lattice evaluation 

The nucleon interpolating field in the physical basis J{x) = e"^'^ \u""^ {x)"^ y$d^ {x)\ u'^{x) can 
be written in the twisted basis at maximal twist as J{x) = -^[^ + iy5\e"^'^ \u'^^ {x)'i§'ysd'^ {x)\ u'^{x). 
The transformation of the electromagnetic current, V^{x), to the twisted basis leaves the form of 
Vii''^ (x) unchanged. We use the Noether lattice current and therefore the renormalization constant 
Zy = 1. The axial current also has the same form in the two bases. In this case we use the 
local current and therefore we need the renormalization constant Za- The value of Za = 0.76(1) [^] 
was determined non-perturbatively in the RF-MOM scheme. This value is consistent with a recent 
analysis [^, which uses a perturbative subtraction of G((p-) terms [|l^] for a better identification 
of the plateau yielding a value of Za = 0.768(3) [^. In order to increase overlap with the proton 
state and decrease overlap with excited states we use Gaussian smeared quark fields | [TT| ] for the 
construction of the interpolating fields: q''{t,x) = 'LyF"^{x,y\U{t)) q^{t,y) with F = {1 + aH)" 
and H{x,y;U{t)) =Y4=iWi{x)5xy-T + Uj{x — i)5x,y+j]. In addition we apply APE-smearing to the 
gauge fields entering H. The smearing is the same as for our calculation of baryon masses with 
the smearing parameters a and n optimized for the nucleon ground state [^]. 

To set the scale we use the nucleon mass in the physical limit. We show in Fig. 2 results at 
three values of the lattice spacings corresponding to j8 = 3.9, j8 = 4.05 and j8 = 4.2. As can be 
seen, cut-off effects are negligible and we can therefore use continuum chiral perturbation theory 
to extrapolate to the physical point. We correct for volume dependence coming from pions propa- 
gating around the lattice [|l^]. To chirally extrapolate we use the well-established &{p^) result of 
heavy baryon chiral perturbation theory (HB;^PT) given by 

3^2 

mN = m%-4cim\--—^m\. (2.1) 

We perform a fit to the volume corrected results at j3 = 3.9, jS = 4.05 and /3 = 4.2 and extract 
ro = 0.462(5) fm. Fitting instead to the j8 = 3.9, jS = 4.05 results we find ro = 0.465(6) fm show- 
ing that indeed cut-off effects are small. To estimate the error due to the chiral extrapolation we use 
HB;^PT to iff{p'^), which leads to ro = 0.489(11). We take the difference between the 0'{p^) and 
^{p^) mean values as an estimate of the uncertainty due to the chiral extrapolation. Fits to other 
higher order ;^PT formulae shown in Fig.3 and described in Ref. [Q] are consistent with &{p^) 
HB;^PT. Using ro = 0.462(5)(27) and the computed ro/a ratios we obtain ap=i, g = 0.089(1)(5), 
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fl'j3=4.05 = 0.070(1)(4) and ap^4 2 = 0.056(2)(3). These values are consistent with the lattice spac- 
ings determined from /;;: and will be used for converting to physical units in what follows. We note 
that results on the nucleon mass using twisted mass fermions agree with those obtained using other 
lattice 0'{a^) formulations for lattice spacings below 0.1 fm ||6|]. 
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Fig. 3: Nucleon mass in units ro. The solid and 
dashed lines are fits to 6{p^) and ff^p^) HB;t;PT. 



Fig. 2: Nucleon mass in units of ro at 3 lattice spac- 
ings and spatial lattice size L such that m-^L > 3.5. 

In order to calculate the aforementioned nucleon matrix elements we calculate respectively the 
two-point and three-point functions: G{q,tf) = Y,^^ e^'^fi T^^ (7« {tf,Xf)lp (0)) and (F^ ,q,t) = 



Y,x,xf ^'''''■'^r^o; (■^«(V'-^/)/(^-^)-^j3(0))> where the projection matrices = 5(1 + 70) and r'= = 
iY^YsYk- The kinematical setup that we used is illustrated in Fig. 1: We create the nucleon at tt = 0, 
f = (source) and annihilate it at tf/a = 12, Pf = (sink). We checked that the sink-source time 
separation of 12a is sufficient for the isolation of the nucleon ground state by comparing the results 
with those obtained when tf /a = \4 is used [0]. We insert the current at t carrying momentum 
q = —pj. In this work we limit ourselves to the calculation of the connected diagram which in the 
isospin limit yields the isovector electromagnetic form factors. This is calculated by performing 
sequential inversions through the sink so that no new inversions are needed for different operator 
j^{t,q). However new inversions are necessary for a different choice of the projection matrices 
r". In this work, we consider the four choices given above, which are optimal for the form factors 
considered here and construct the ratio 



^ G^(r,g,0 G{pi,tf-t)G{0,t)G{0,tf) tf^u^^^^ 

GiO,tf) \ G{0,tf-t)Gipi,t)G{pi,tf) ~^ ^ '^^ ■ 



(2.2) 



The leading time dependence and overlap factors cancel yielding as the plateau value n'^(r, 
from which we extract the form factors using the relations 

n^r",^) = ^[{m+E)do,^ + iqtdt.^] Ge{q^), n'(r\^) = ^^Z^juqjSijGMiQ^) 



2m 



J' 



2nP- 
E(E+m) 



and U^i{T\q) = ^^[^Gp{Q^)-{E + m)8aGA{Q^)\, k=\,---,3, where c = 
3. Results 

3.1 Isovector electromagnetic form factors 

The isovector electric and magnetic form factors at jS = 3.9 are shown in Figs. 4 and 5 respec- 
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lively on lattices of size 24^ x 48 and 32^ x 64 and for the pion masses and statistics listed in Fig. 6. 
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Fig. 4: The isovector electric form factor as a func- 
tion of 2^ compared to experiment (solid curve). 



(GeV^) 

Fig. 5:The isovector magnetic form factor as a func- 
tion of compared to experiment (solid curve). 
By fitting the isovector magnetic form factor to a dipole form we extract the isovector magnetic 
moment and Dirac and Pauli root mean squared (r.m.s.) radii shown in Figs. 6 and 7. As can be 
seen, the results obtained using twisted mass fermions are in agreement with recent results using 



dynamical domain wall fermions (DWF). Heavy baryon chiral perturbation theory to one-loop [ |13| ] 
can be used to extrapolate to the physical point. In the case of the Dirac r.m.s radius we fit the 
product r^Ky so that only one fit parameter enters just as in the case of r\. This shifts the curve but 
does not affect its slope. We show fits to our results alone as well as when we include the results 
obtained by the RBC-UKQCD collaborations [Q]. The magnetic moment with three fit parameters 
is reproduced whereas for r\ we obtained a weaker dependence on the pion mass as compared to 
the one predicted in chiral expansions. 
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Fig 7: The Dirac (upper) and Pauli (lower) r.m.s 
radii. 
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3.2 Axial charge 

Our results on the nucleon axial charge are shown in Fig. 8 and are in agreement with those 
obtained using domain wall fermions. Within our errors no sizable finite volume effects are ob- 
served. In order to extrapolate to the physical point we use one-loop chiral perturbation theory in 
the small scale expansion 



There are three parameters to fit: 
gA(0), the value of the axial charge 
at the chiral point, g\ and a counter- 
term CssE- Fitting using the TMF re- 
sults we findgA(O) = 1.10(13), gi = 
6.05(4.66) wACssE = -3.65(3.58). 
The parameters g\ and Csse are 
highly correlated explaining the re- 
sulting large error band. Including in 
the fit the DWF data that span larger 
pion masses we obtain a very differ- 
ent curve, showing the sensitivity in 
the chiral extrapolation of gA- 

3.3 Axial form factors 
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Fig. 8: The nucleon axial charge. The solid (dashed) curve is a 
chiral fit using TMF (TMF and DWF) results. 



Results for the axial form factors Ga{Q^) and Gp{Q^) are shown in Figs. 9 and 10 respectively. 
We perform a dipole fit to Ga(G^) extracting an axial mass larger than in experiment as expected 
from the smaller slope shown by the lattice data both for TMF and DWF. Assuming pion pole 
dominance we can relate the form factor Gp{q^) to Ga{Q^)- Using the pion mass measured on the 
lattice we predict the dashed curve shown in Fig. 10. Our lattice data on Gp{q^) are flatter than 
pion pole dominance predicts requiring a larger pole mass than the pion mass. Large volume effects 
are expected at low indicated by the deviation of Gp{Q^) from the fitted curve at the smallest 
2^-value. 

20 
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Q2 (GeV^) 

Fig. 9: Axial form factor Ga[Q^) as a function of 
Q^. The dashed (solid) line is the best dipole fit to 
the lattice (experimental) results. 
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Fig. 10: The dotted line is a fit of lattice results to 
the form ^(^q1^„i-^ and the dashed (solid) line is the 
prediction of pion pole dominance using lattice (ex- 
perimental) results on Ga{Q^)- 
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4. Conclusions 

Using Nf = 2 twisted mass fermions we obtain accurate results on the isovector electromag- 
netic Ge,Gm and axial GA,Gp form factors as a function of for pion mass in the range of about 
260-470 MeV. The general feature is a flatter dependence on than experiment. The Dirac r.m.s 
radius thus shows a weaker dependence on the pion mass than expected from chiral perturbation 
theory. Finite volume effects are found to be small on quantities like gA and the isovector magnetic 
moment and r.m.s radii. Our results are in agreement with recent results obtained using dynamical 
A^F = 2 + 1 DWF. At the physical point using TMF we find = 1 • 1 3 ( 1 0) close to the experimental 
value albeit with a large error due to the chiral extrapolation. An analysis of these form factors at 
J3 = 4.05 is under way so that a check of cut-off effects can be carried out. 
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